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Understanding the DNA elements that constitute and control the regulatory genome is critical for the appropriate 
therapeutic management of complex diseases. Here, using chromosome Y (ChrY) consomic mouse strains on the C57BL/ 
6] (B6) bacicground, we show that susceptibility to two diverse animal models of autoimmune disease, experimental 
allergic encephalomyelitis [EAE] and experimental myocarditis, correlates with the natural variation in copy number of 
Sly and Rbmy multicopy ChrY genes. On the B6 background, ChrY possesses gene regulatory properties that impact 
genome-wide gene expression in pathogenic CD4^ T cells. Using a ChrY consomic strain on the S]L background, we 
discovered a preference for ChrY-mediated gene regulation in macrophages, the immune cell subset underlying the EAE 
sexual dimorphism in S]L mice, rather than CDV T cells, importantly, in both genetic backgrounds, an inverse correlation 
exists between the number of S/y and Rbmy Chr\ gene copies and the number of significantly up-regulated genes in immune 
cells, thereby supporting a link between copy number variation of S/y and Rbmy w/\th the ChrY genetic element exerting 
regulatory properties. Additionally, we show that ChrY polymorphism can determine the sexual dimorphism in EAE and 
myocarditis, in humans, an analysis of the CD4^ T cell transcriptome from male multiple sclerosis patients versus healthy 
controls provides further evidence for an evolutionarily conserved mechanism of gene regulation by ChrY. Thus, as in 
Drosophila, these data establish the mammalian ChrY as a member of the regulatory genome due to its ability to epige- 
netically regulate genome-wide gene expression in immune cells. 



[Supplemental material is available for this article.] 

Sex-specific differences exist in many aspects of immune system 
physiology and contribute to the pathogenic differences in in- 
fectious disease observed between males and females (Pennell et al. 
2012). Females tend to mount a more robust immune response 
against many infectious pathogens, leading to a better prognosis in 
disease outcome compared with males. However, the evolutionary 
advantage of this heightened immune response in females also 
contributes to their higher risk in developing autoimmune (AI) 
disease. While these sex differences in immunity are predominately 
linked to the differential effects of sex hormones on immune cells, 
genetic differences in sex chromosomes can also influence the im- 
mune response and susceptibility to disease (Ober et al. 2008; 
Arnold and Lusis 2012). Importantly, although the sex chromo- 
some effect is often attributed to the difference in the number of 
X chromosomes (ChrX) between males and females, increasing 
evidence indicates that the mammalian Y chromosome (ChrY) 
contributes to immune responses and susceptibility to AI and in- 
fectious diseases (Teuscher et al. 2006; Spach et al. 2009; Robinson 
et al. 2011; Case et al. 2012; Charchar et al. 2012). 

Like autosomal chromosomes, ChrX and ChrY are thought to 
have once been identical pairs that were free to recombine and 
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exchange genetic material. Over the course of evolution, ChrY 
became unique from all other chromosomes due to the acquisition 
of a dominant sex-determining gene and subsequent chromo- 
somal inversions that restricted recombination and led to its deg- 
radation (Charlesworth 1991; Gschwend et al. 2012). Although 
ChrY has evolved independently across different taxa, they share 
common genetic features that are specific to ChrY. For example, 
the relatively few protein-coding genes ChrY possesses are pre- 
dominately male-specific genes acquired through transposition 
and translocation from other chromosomes (Skaletsky et al. 2003; 
Koerich et al. 2008). The remainder of ChrY is largely hetero- 
chromatic and composed of multicopy genes, repeat sequences, 
and transposable elements. Therefore, the nature of ChrY has led 
to the consensus that it is primarily composed of "junk" DNA 
whose contribution to phenotypic differences among the sexes is 
limited to sexual development and spermatogenesis. Indeed, until 
recently, large deletions in the long arm of the murine ChrY had no 
overt phenotypic effects except on spermatogenesis and sex ratio 
distortion (Tiepolo and Zuffardi 1976; Styrna et al. 1991a,b; Sun 
et al. 2013). 

The first phenotypic evidence describing alternative func- 
tions for ChrY emerged from Drosophila studies linking ChrY 
polymorphism with behavior phenotypes (Stoltenberg and Hirsch 
1997; Huttunen and Aspi 2003). Aside from its role in male de- 
velopment, the Drosophila ChrY differentially regulates the ex- 
pression of autosomal and ChrX genes, including a preferential 
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enrichment of genes associated with the immune response, by 
altering chromatin between the heterochromatic and euchro- 
matic state, thereby influencing numerous complex phenotypes 
(Lemos et al. 2008, 2010; Piergentili 2010; Paredes et al. 2011). 
The ability of ChrY to alter chromatin dynamics and regulate 
gene expression is attributed to natural polymorphic variation in 
the multicopy ribosomal genes present on ChrY (Paredes et al. 
2011; Zhou et al. 2012). 

The mammalian ChrY can also moderate complex pheno- 
types not typically associated with this chromosome (Monahan 
and Maxson 1998; Chen et al. 2013). In mice, genetic variation in 
ChrY influences susceptibility of male C57BL/6J (B6) mice to ex- 
perimental allergic encephalomyelitis (EAE), the animal model of 
multiple sclerosis (MS), as well as the EAE sexual dimorphism ob- 
served in the SJL/J mouse strain (Teuscher et al. 2006; Spach et al. 
2009). Natural variation in ChrY also impacts the mortality rate of 
B6 mice following infection with coxsackievirus B3 (CVB3) and 
the basal levels of invariant natural killer T cells independently of 
the sex determining region of ChrY (Sry) expression level poly- 
morphisms or differences in testosterone levels (Case et al. 2012). 
Furthermore, in humans, a recent study on coronary artery disease 
revealed a link between disease susceptibility and Haplogroup I of 
the male-specific region of human ChrY (Charchar et al. 2012). 
However, the mechanism whereby natural genetic variation in the 
mammalian ChrY influences susceptibility to disease is unknown 
and is thus the focus of this study. 

Here, we show that susceptibility to two diverse animal 
models of Al disease correlates with the natural variation in copy 
number of the ChrY multicopy genes Sycp3 like Y-linked (Sly) and 
RNA-binding motif, ChrY (Rbmy). ChrY possesses gene regulatory 
properties that impact both genome-wide gene expression and 
alternative splicing in a cell type-specific manner depending on 
the genetic background of the mice. Furthermore, an inverse cor- 
relation exists in the number of Sly and Rbmy ChrY multicopy 
genes and the number of significantly up-regulated genes in im- 
mune cells, thereby supporting a link between copy number var- 
iation of Sly and Rbmy with the ChrY genetic element exerting 
regulatory properties. Additionally, by comparing female B6 mice 
with B6-ChrY consomic lines we show that ChrY polymorphism 
also determines the gender specificity of the sexual dimorphism 
in Al disease. Thus, these data support an evolutionarily con- 
served mechanism retained by mammalian ChrYs to influence 
complex phenotypes, including susceptibility to Al disease, by act- 
ing as a global trans-expression quantitative trait locus (trans-eQTL) 
in males. 

Results 

ChrY consomic strains exhibit a continuous distribution 
in Al disease severity characteristic of a quantitative trait 

Investigating the genetic components underlying ChrY's ability to 
regulate susceptibility to Al disease in male mice is complicated by 
its structure, which prevents traditional methods of candidate 
gene identification. However, linkage of disease susceptibility to 
specific ChrY polymorphisms can be achieved by comparing the 
strain distribution pattern (SDP) of disease susceptibility with the 
SDPs for natural genetic variations segregating among ChrYs. This 
comparison is made possible through the use of ChrY consomic 
strains of mice, also known as chromosome substitution strains. 
These mice are generated by intercrossing a male mouse possess- 
ing a ChrY of interest to a female mouse with the desired genetic 



background, which in this study is B6. Then, through a series of 
backcrosses of 10 generations or more to B6 females, the auto- 
somes, ChrX, the pseudoautosomal (PAR) regions of ChrX and 
ChrY, and the mitochondrial genome of the strain donating ChrY 
are replaced with the B6 genome. Therefore, the genetic variation 
among ChrY consomic strains presumably comes from the non- 
recombining region of ChrY called the nonpseudoautosomal or 
male-specific region (Nadeau et al. 2012). 

To identify the source of the ChrY polymorphism controlling 
susceptibility to Al disease, we monitored the severity of both EAE 
and CVB3-induced myocarditis among a panel of B6-ChrY con- 
somic strains, where the mouse strain donating ChrY to B6 is in- 
dicated in superscript (e.g., B6-ChrY^^J). To test for susceptibility 
to EAE, adult male mice from 10 different B6-ChrY consomic 
strains were immunized with 100 [xg of myelin oligodendrocyte 
glycoprotein peptides 35-55 (MOGss.ss) in complete Freund's 
adjuvant (CFA) on day 0 and day 7 (2x protocol) and scored for 
clinical signs of EAE for 30 d following immunization. Analysis of 
the severity of the clinical disease course (P < 0.00001) (Fig. lA) 
and the cumulative disease scores (CDS) (overall P < 0.0009) 
(Fig. IB) for B6 and each of the B6-ChrY consomic strains re- 
vealed a continuous distribution in disease severity across the 
strains, consistent with quantitative inheritance (Castle 1914; 
East 1916). The severity of the clinical disease course for each 
consomic strain was compared with B6 identifying B6-ChrY^^J, 
B6-ChrY^^^, B6-ChrY™^, B6-ChrY^^, B6-ChrY^u^, and B6- 
ChrY^^^ as being equivalently susceptible to EAE as B6 (Supple- 
mental Fig. SIA). In contrast, B6-ChrY^^ B6-ChrY^^^, B6-ChrY^^, 
B6-ChrY^^^^^ and B6-ChrY^J^ showed a significant reduction in 
the severity of the course of clinical disease compared with B6 
(Supplemental Fig. SIB). Among these low-responder strains, B6- 
ChrY^J^ exhibited the largest difference in both clinical disease 
and neuropathology (Fig. IC). 

Myocarditis is an inflammatory heart disease with an Al 
component that is often caused by local infections and pre- 
dominates in both men and male experimental mice (Gauntt and 
Huber 2003; Fairweather et al. 2013). We recently reported an ef- 
fect of sex chromosome complement on coxsackievirus B3 (CVB3)- 
induced myocarditis using the 'Tour core genotypes" (FCG) mouse 
model designed to independently investigate the contribution of 
sex chromosomes (XX vs. XY) or gonadal type (ovaries vs. testes) 
on disease (Arnold and Chen 2009). Compared with gonadecto- 
mized XX-female and XX-male mice, gonadectomized XY-female 
and XY-male mice exhibited less severe disease, indicating that 
myocarditis susceptibility is influenced by ChrX and/or ChrY 
(Robinson et al. 201 1). The response of FCG mice to infection with 
CVB3 and the development of CVB3 -induced myocarditis in XX 
animals are similar to the response of SJL-FCG mice to EAE, in 
which the XX sex chromosome complement is associated with 
significantly more severe disease than the XY complement (Smith- 
Bouvier et al. 2008). The interpretation of this outcome favored by 
these authors is that the increase in EAE severity seen in female SJL 
mice is primarily a function of the XX complement. However, 
using SJL-ChrY^^^'^ consomic mice, we also show that ChrY^^^ can 
influence the difference in EAE susceptibility between female 
(X^J^X^J^) and male (X^J^Y^J^) mice (Spach et al. 2009). Therefore, to 
assess whether susceptibility to myocarditis is similarly regulated 
by natural polymorphic variation in ChrY, male mice from four 
B6-ChrY consomic strains were infected i.p. with 50 PFU CVB3, 
and their hearts were evaluated for inflammation and cardiac virus 
titers 7 d after infection. We found that B6-ChrY^JS B6-ChrY^^ 
and B6-ChrY^™ mice, but not B6-ChrY^^^ mice, exhibit a signif- 
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icant decrease in CVB3-induced myocarditis relative to B6 (Fig. 
ID). Residual heart virus titers did not predict the heart score 
among the strains with the exception of B6-ChrY^^^, which had 
both significantly lower viral titers (P < 0.05) as well as lower 
myocarditis scores (P < 0.001) compared with B6 (Fig. 1D,E and 
data not shown). Taken together, two diverse models of AI disease 
exhibit a continuous distribution in disease severity across the 
B6-ChrY consomic strains, which is characteristic of quantitative 
inheritance. 
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Sry expression level polymorphisms and variation 

in testosterone production do not correlate with EAE severity 

Previously we showed that neither the serum testosterone levels 
nor the survival of CVB3-infected B6-ChrY male mice were con- 
trolled by Sry expression level polymorphisms (Case et al. 2012), 
which when inherited by B6 mice heterozygous for the t-complex 
encoded TOrleans mutation on Chrl7, leads to varying degrees of 
sex reversal (Nagamine et al. 1999; Washburn et al. 2001; Albrecht 
et al. 2003). Similarly by comparing the EAE-CDS for B6 and each 
consomic strain with their serum testos- 
terone levels, we found that EAE severity 
also does not correlate with differences in 
testosterone production (r = 0.36; P = 
0.27), a hormone that is known to in- 
fluence EAE and myocarditis severity in 
males (Fig. IF; Huber et al. 1982; Lyden 
et al. 1987; Bebo et al. 1998; Fillmore et al. 
2004; Frisancho-Kiss et al. 2009; Spence 
andVoskuhl 2012). 

To examine whether Sry expression 
level polymorphisms affect susceptibility 
to AI disease, we compared the SDP of Sry 
alleles with EAE-CDS across the B6-ChrY 
consomic strains. This revealed that Sry 
polymorphisms similarly do not influ- 
ence disease severity, suggesting that Sry 
expression polymorphisms do not under- 
lie the differences in disease phenotypes 
(P = 0.165) (Fig. IG). Furthermore, as pre- 
viously reported, we found that Sry poly- 
morphisms do not influence the level of 
adult testosterone production in the B6- 
ChrY consomic strains (P = 0.354) (Case 
et al. 2012). Additionally, an analysis of 
the single nucleotide polymorphism (SNP) 
variation using the Mouse Phenome 
Database identified only two ChrY haplo- 
groups among the strains with limited 
diversity that are also discordant with 
disease susceptibility (Supplemental Table 
S2). Given the limited SNP diversity be- 
tween the annotated regions on the short 
arm of ChrYs, it is important to consider 
the impact that other polymorphisms, 
such as structural polymorphism, may 
have on regulating phenotypes. 
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ChrY gene copy number variation 
correlates with EAE and myocarditis 
severity among the B6-ChrY consomic 
strains 

Structural polymorphism can arise through 
variations in the number of repeat se- 
quences, inverted sequences, and ret- 
roelements present in ChrY DNA. In 
Drosophila, ChrY-mediated phenotypic 
differences are attributed to natural poly- 
morphic variation in multicopy ribosomal 
genes present on ChrY (Lemos et al. 2008, 
2010; Paredes et al. 2011; Zhou et al. 2012). 
The musculus ChrY possesses several 
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multicopy genes, including Sly, Rbmy, and spermiogenesis tran- 
script on the Y 1 (Sstyl), that differ in number between ChrYs of 
inbred strains (Fig. 2A). Therefore, to investigate whether copy 
number variation in these multicopy ChrY genes may impact 
susceptibility to EAE, we ran a correlation analysis between gene 
copy number and CDS. This analysis shows that the number of 
gene copies of Sly and Rbmy, but not Sstyl, correlates with the 
severity of EAE (Fig. 2B). Similarly, a correlation analysis between 
gene copy number and myocarditis heart score also reveals a striking 
association between the number of gene copies of Sly and Rbmy, but 
not Sstyl, with the severity of CVB3 -induced myocarditis (Fig. 2C). 

The mammalian ChrY acts as a trans-eQTl shaping immune cell 
transcriptomes 

Given that the ChrY-mediated influence on AI disease suscepti- 
bility is indicative of a quantitative trait, we sought to determine 
whether ChrY can act as frans-eQTL and regulate global gene 
expression in cells central to disease pathogenesis (Bhasin et al. 
2008). We analyzed the transcriptome of naive CD4^ T cells, which 
are essential to disease pathogenesis in B6 mice, between resistant 
B6-ChrY^^^ and susceptible B6 male mice by microarray analysis. 
RNA was isolated from fluorescently activated cell sorted (FACS) 
CD4"^TCRp^ cells obtained from the lymph nodes of adult (8 wk of 
age) naive B6-ChrY^J^ and B6 male mice, and mRNA expression, the 
presence of alternative splice variants, and miRNA expression were 
assessed in each sample. 

After filtering the gene array data using a false discovery rate 
(FDR) <0.05, we identified 734 differentially expressed transcripts 
between B6-ChrY^J^ and B6 CD4+ T cells. Furthermore, 3247 tran- 
scripts were identified as being alternatively spliced, with 398 of 
these being both alternatively spliced and differentially expressed 
(Fig. 3A). These differentially regulated genes are distributed 
throughout the genome and reside on all autosomes and both sex 
chromosomes. Interestingly, no miRNA probe sets passed an FDR 
of 0.05. These data demonstrate the magnitude of ChrY's ability to 
act as a frans-eQTL and epigenetically regulate the transcriptome, 
particularly in relation to alternative splicing, and raise the possi- 



Figure 1. Natural genetic variation in ChrY mediates susceptibility to autoimmune disease. (A) Male 
mice from the B6-ChrY consomic strains were immunized with MOG35_55 using the 2x protocol and 
the clinical score was monitored over 30 d. The consomic strain represented by each line is color-coded 
with the bar graph in B. Wild-type B6 male are represented by black circles. The significance of the 
differences in disease course among the strains was determined by two-way ANOVA (interaction [F = 
1.18; DFn = 319; DFd = 7770; P = 0.02], day post-injection [f = 132.81; DFn = 29; DFd = 7770; P< 
0.0001], and strain [F= 23.53; DFn = 1 1; DFd = 7770; P< 0.0001]). (B) The degree of EAE susceptibility 
among the B6-ChrY consomic mice was calculated by subtracting the clinical disease score (CDS) of B6 
from the CDS of each strain. The x-axis indicates the strain origin of the Y chromosome. F= 2.58, P < 
0.0009. A heterogeneity test was used between cohorts of mice with no significant differences detected. 
Thus, data were pooled from three independent experiments and the total number of animals analyzed 
in A and B are indicated in Supplemental Table SI . (C) Immunohistopathology of the brain and spinal 
cord from B6-ChrY consomic strains immunized using the 2x protocol. Only the total histopathology 
score is shown. Significance determined by one-way ANOVA. (**) P < 0.01 . (D) Male mice from wild- 
type B6 and four consomic lines were infected with 50 PFU CVB3 and their hearts were evaluated for 
myocarditis. Heart score of B6 was subtracted from the heart score for each consomic. Labeling on the 
X-axis represents the strain donating ChrY to B6. Results are representative of three individual experi- 
ments, n > 5 mice per strain. Significance determined by one-way ANOVA and Dunnett's multiple 
comparison test. (*) P < 0.05; (**) P < 0.01 ; (***) P < 0.001 . (£) Linear regression analysis of heart score 
and heart viral titers reveals the absence of a significant relationship between the regressed variables. Left 
X-axis (orange linear regression line and orange circles) corresponds to the B6-ChrY consomic data 
combined from each of the four tested strains, f =0.72, P=0.41 . Right x-axh (black linear regression line 
and black squares) corresponds to the WT B6 data, f = 0.43, P = 0.52. (F) Correlation analysis of serum 
testosterone with the CDS for EAE (r= 0.36; P= 0.27). (C) Functionally significant 5Ay polymorphisms are 
discordant with ChrY-mediated differences in disease susceptibility. Data from B were grouped 
according to their functional Sry polymorphisms. 



bility that gene regulation may be independent of miRNA expres- 
sion in CD4"'T cells. 



Differential expression of chromatin remodeling genes 
in B6-ChrY^^'- vs. B6 CD4^ T cells 

We utilized the Mouse Genome Informatics (MGI) Gene Ontology 
(GO) Term Finder to identify the GO terms containing a significant 
enrichment of genes within the B6-ChrY^J^ vs. B6 alternatively 
spliced and differentially expressed transcript lists (Fig. 3A; Sup- 
plemental Tables S3, S4). Both lists contained multiple terms as- 
sociated with chromatin dynamics whose dysregulation has the 
potential to lead to the ChrY-dependent genome-wide transcrip- 
tome differences observed in CD4"^ T cells. Therefore, to confirm 
these gene associations, we used a quantitative PGR (qPCR) array 
specific for mouse epigenetic chromatin remodeling factors to as- 
sess their differential expression in B6-ChrY^J^ vs. B6 CD4^ T cells. 
Of the 84 chromatin remodeling genes represented on the qPCR 
array, 54 exhibited significant differences in gene expression 
between B6-ChrY^J^ and B6 004"" T cells, consistent with ChrY- 
mediated changes in chromatin remodeling gene expression 
identified in the microarray (Fig. 3B; Supplemental Table S5). 

ChrY-mediated transcriptome differences translate to 
alterations in cytokine production in stimulated CD4"^ T cells 

We used Ingenuity Pathway Analysis (IPA) to identify specific ca- 
nonical pathways that might be affected by the ChrY-dependent 
changes in the CD4^ T-cell transcriptome. We included the genes 
from both the alternatively spliced and differentially expressed 
lists and used the core analysis feature of IPA. The top 10 canonical 
pathways identified by IPA include those involved in gene ex- 
pression, protein synthesis, T- and B-cell signaling, and glucocor- 
ticoid receptor signaling, consistent with the MGI GO term data 
(Supplemental Fig. S3). Furthermore, three of the pathways were 
specific for signaling in T cells, including T-cell receptor signaling, 
leading us to hypothesize that the activation profile should 
markedly differ between B6-ChrY^J^ and B6 CD4+ T cells (Fig. 3C). 

To address this, purified CD4^ T cells from 
B6-ChrY^J^ and B6 lymph nodes were 

stimulated with anti-CD3 and anti-CD28 

monoclonal antibodies and the superna- 
tants were tested for IL2, IFNG, IL17A, 
TNF, and IL6 production after 24, 48, and 
72 h by ELISA. We observed significant 
reductions in IFNG, IL17A, and IL6 cy- 
tokine secretion by stimulated CD4"^ 
T cells from B6-ChrY^J^ compared with 
B6 (Fig. 3D). These data indicate that the 
ChrY-mediated changes in the CD4"^ T-cell 
transcriptome, as identified by microarray 
analyses, translates directly into altered 
protein synthesis as a result of biological 
changes in T-cell activation and signaling. 



ChrY exerts cell type-specific effects 
on gene regulation depending on 
the autosomal background of the mice 

Unlike the B6 strain, SJL mice exhibit 
a sexual dimorphism in EAE that is the 
result of an age-dependent resistance in 
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Figure 2. 5// and Rbmy gene copy number correlates with EAE and CVB3 susceptibility among B6-ChrY consomic strains. (A) Gene copy number for 
each of the consomic strains relative to B6. n > 5 mice per strain. (B) Correlation analysis of gene copy number with the CDS for EAE. (C) Correlation 
analysis of gene copy number with myocarditis heart score. The x-axis represents the strain origin of ChrY. 



males due, in part, to the inability of macrophages from young 
SJL mice to efficiently process and present antigen to effector cells 
(Stohlman et al. 1985; Cua et al. 1995; Papenfuss et al. 2004). 
Furthermore, this phenotype is reliant on the ChrY^^^, because 
young SJL male mice inheriting ChrY^^°-^ (SJL-ChrY^^^-^) are as 
susceptible to disease as older males (Spach et al. 2009). To test 
whether ChrY acts as a frans-eQTL in the immune cells essential to 
disease susceptibility in SJL mice, we analyzed the transcriptomes 
of FACS-purified CD4"^ T cells and adherent macrophages from 
young (<4 wk of age) and old (>6 mo of age) male SJL-ChrY^^°"^ 
and SJL mice. An analysis of the main effect of strain on each of 
the cell types suggest that ChrY has a greater influence on the 
transcriptome of macrophages compared with CD4^ T cells, with 
2932 vs. 42 differentially expressed transcripts passing an FDR < 
0.05, respectively (Fig. 4A,B). 

Similar to the B6-ChrY consomic comparison, ChrY also 
greatly influences the alternative splicing of genes in the SJL- 
ChrY^^°-^ and SJL strain combination, with 2987 vs. 348 alter- 
natively spliced transcripts in macrophages and CD4"^ T cells, 
respectively (Fig. 4A,B). MGI GO term analysis on the CD4^ T-cell 
data set indicated that the alternatively spliced transcripts con- 
tained significant term associations for transcriptome regulation 
whereas no significant term enrichments were identified for the 



differentially expressed transcript list (Fig. 4A; Supplemental 
Table S6). Interestingly, the genes represented in the macrophage 
array have very different signatures depending on whether they 
are alternatively spliced or differentially expressed. Alternatively 
spliced genes typically function in regulating gene expression, 
as was observed for CD4^ T cells, whereas the differentially ex- 
pressed genes function in other diverse cellular processes un- 
related to gene regulation (Fig. 4B; Supplemental Tables S7, S8). 
These data reveal that allelic differences in ChrY affects gene 
regulation more in macrophages than CD4^ T cells in the SJL 
genetic background. Thus, these data provide further support for 
the existence of gene regulatory elements within ChrY that exert 
regulatory properties and show that these effects are cell-type- 
and autosomal background-dependent. 

The number of S/y and Rbmy gene copies is inversely correlated 
with the up-regulation of gene expression 

The microarray analyses described above led to a striking observa- 
tion in the directionality of transcript expression within both ChrY 
consomic backgrounds. In B6-ChrY^J^ CD4^ T cells, >75% of the 
differentially expressed transcripts were preferentially up-regulated 
compared with wild-type B6. In contrast, placing ChrY^^°-^ on the 
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Figure 3. ChrY-nnediated regulatory variation in genome-wide gene expression influences the activation profile of CD4^ T cells. (A) Venn diagram 
representing the number of genes passing an FDR < 0.05 between B6-ChrY^''- and B6 CD4^ T cells. MGI GO Term Finder was used to identify enriched 
biological processes associated with each gene list, examples of which are shown. (B) Heat map from quantitative PGR array specific for mouse epigenetic 
chromatin remodeling factors on three biological replicates of purified GD4^ T cells from B6 and B6-GhrY^''-. Red represents up-regulated genes and green 
down-regulated genes. Degree of change indicated by intensity of color. (C) IPA TGR signaling pathway with the molecules differentially regulated 
between B6-GhrY^^'- and B6 GD4^ T cells highlighted in orange. (D) Purified GD4^ T cells from B6-GhrY^''- and B6 lymph nodes were incubated with anti- 
GD3 and anti-GD28 monoclonal antibodies and the supernatants were tested for IL2, IFNG, IL1 7 A, TNF, and IL6 production after 24, 48, and 72 h by ELISA. 
Data are representative of three independent experiments. Significance determined by two-way ANOVA followed by Bonferroni post-hoc multiple 
comparison test. (*) P < 0.05; (****) P < 0.0001 . 



SJL background resulted in the down-regulation of >85% of the 
differentially expressed transcripts in CD4^ T cells and —60% of 
the differentially expressed transcripts in macrophages compared 
with the wild-type ChrY^^^. Therefore, ChrY^^^ contains an overall 
low number of multicopy genes and an increase in the number of 
up-regulated genes. Conversely ChrY^^^-^ contains a high number 
of multicopy genes, which is equivalent to that of ChrY^^, and 
a decrease in the number of up-regulated genes (Fig. 4C). These 
data indicate that an inverse correlation exists between Sly and 
Rbmy copy number variation and the up-regulation of gene ex- 
pression in immune cells, thereby providing a link between the 
copy number variation of Sly and Rbmy with the ChrY frans-eQTL 
regulatory properties. 



Evidence for the evolutionary conservation of ChrY 
as an eQTL in MS susceptibility 

A recent study focusing on coronary artery disease in humans re- 
vealed a link between disease susceptibility and ChrY haplogroups 
in association with distinct autosomal and ChrX transcriptional 
profiles in macrophages (Charchar et al. 2012). This work indicates 
that, like in Drosophila and mice, ChrY polymorphism may in- 
fluence disease susceptibility by regulating the gene expression 
profiles of pathogenic immune cells in men. Therefore, to identify 
pathways known to be directly regulated by ChrY polymorphism 
in mice and that are evolutionarily conserved in humans, we uti- 
lized a naive CD4^ T-cell expression data set from male MS patients 



Genome Research 1479 



www.genome.org 



Case et al. 



CD4*T Cells 

SJL-ChrY"<^^vsSJL 
Alternaf vely Spliced Differentially expressed 



Macrophages 

SJL-ChrY^i°^ ^5 SJL 




Alternatively Spliced 



Differentially expressed 




Gene Expression 
Chromatin Modification 
RNA Processing 



Gene Expression 
Chromatin Modification 
RNA Processing 



Response to Stress 
Inflammatory Response 
Cell Chemotaxis 



Upregulated Genes 
CNV 



CD4* T cells 

Male CIS BG-ChrYS^^ 



Healthy Male Controls 




SJL(CD4+) BIO.S(Macs) 
Strain Origin of ChrY 




Gene Expression, Transcription, Chromatin 
fy/lodification, Immune Response 



Figure 4. ChrY preferentially regulates the transcriptome of macrophages compared with CD4^ T 
cells in the SJL genetic background. Venn diagram representing the number of genes passing an FDR < 
0.05 between SJL-ChrY^^°-^and SJL. (A) CD4^T cells and (8) macrophages. MGI GO Term Finder was 
used to identify enriched biological processes associated with each gene list, examples of which are 
shown. (C) Graph representing the percentage of up-regulated genes in the ChrY consomic strain vs. 
the parental in relation to the averaged copy number variation for 5// and Rbmy. (D) A Venn diagram 
representing the genes differentially regulated between the human and B6 mouse CD4^ T-cell micro- 
arrays, 440 of which are shared between the two species. MGI GO term analysis on these shared terms 
identified an enrichment of genes involved in gene regulatory processes. 



sex chromosome complement (Robinson 
et al. 2011). In contrast, B6 mice do not 
typically exhibit a sexual dimorphism in 
EAE (Papenfuss et al. 2004), which is 
consistent with the lack of a sex chro- 
mosome complement effect observed in 
the B6-FCG model (Smith-Bouvier et al. 
2008). 

Interestingly, by comparing the 
myocarditis heart scores of B6-ChrY 
consomic strains with the heart score of 
female B6 mice, the sexual dimorphism is 
eliminated by ChrY^^^ (Fig. 5A). Fur- 
thermore, by comparing the EAE clinical 
disease course of B6-ChrY consomic 
strains with the disease course of female 
B6 mice, we find that varying degrees of 
sexual dimorphism in EAE severity emerge 
in this otherwise non-sexually dimorphic 
strain (Fig. 5B,C). No sexual dimorphism 
in EAE severity was observed between B6 
females and male B6, B6-ChY^/J, B6- 
ChrY^^^, B6-ChrY^^'^, B6-ChrY^^, or B6- 
ChrY^^^ consomic mice. In contrast, 
compared with B6 females, B6-ChrY^^^ 
males exhibited increased disease severity 
while male Be-ChrY^^^^^ B6-ChrY^^, 
B6-ChrY^^, B6-ChrY^JS and Be-ChrY^^'^ 
consomic mice presented with less severe 
disease (Fig. 5C). Clearly, the inheritance 
of polymorphic ChrYs has a profound ef- 
fect on disease susceptibility in B6 males 
and may play an important role in de- 
termining the gender bias observed in the 
Al diseases studied. 



exhibiting the earliest clinical signs of disease, known as clinically 
isolated syndrome (CIS), with healthy controls available through 
the Gene Expression Omnibus (Series GSE13732) (Corvol et al. 
2008). After filtering the data set using a FDR < 0.05, we identified 
1556 differentially expressed genes between male CIS and con- 
trols. Overlapping this filtered data set with the B6-ChrY^J^ vs. B6 
mouse CD4"^ T-cell data set identified 440 genes common between 
mouse and human. This list was specifically enriched for genes 
functioning in gene expression (P= 3.8 X 10~^^), transcription (P = 
6.62 X 10"^°), chromatin modification (P = 4.02 X 10"^^), and 
immune system processes (P = 2.07 X 10~^^) (Fig. 4D; Supple- 
mental Table S9). Therefore, this analysis provides preliminary 
evidence for the evolutionary conservation of ChrY as an eQTL 
capable of modifying the regulatory genome of CD4"^ T cells and 
MS susceptibility in men. 

Genetic variation in ChrY modifies the sexual dimorphism 
in myocarditis and leads to a sex bias in EAE severity 

Autoimmune diseases often exhibit a sexual dimorphism in their 
incidence, clinical course, and/or severity (McCombe et al. 2009). 
MS predominates in women by a 3:1 ratio, whereas myocarditis 
predominates in men by a 2:1 ratio, and in both Al diseases, men 
often develop more severe disease compared with females. Simi- 
larly, CVB3-induced myocarditis is more severe in male B6 mice 
than in females, and this sexual dimorphism is influenced by the 



Discussion 

It is well established that ChrY protein-coding genes are critical for 
the sexual development and fertility of males. Beyond this, how- 
ever, there is only a small, though growing, body of research ex- 
ploring the existence of alternative functions for ChrY in de- 
termining organismal phenotypes (Monahan and Maxson 1998; 
Burgoyne et al. 2002; Wesley et al. 2007; Spach et al. 2009; Arnold 
2012; Case et al. 2012; Charchar et al. 2012; Chen et al. 2013). Here 
we show that natural polymorphic variation in ChrY impacts 
susceptibility to two diverse animal models of Al disease by acting 
as a trans-eQTL in immune cells linked to disease pathogenesis. The 
genome-wide regulatory properties of ChrY, which impact both 
gene expression and alternative splicing, establish this chromosome 
as a member of the regulatory genome in males. Importantly, an 
inverse correlation exists between the copy number variation in Sly 
and Rmby ChrY multicopy genes and the number of significantly 
up-regulated genes in immune cells, thereby providing a link 
between copy number variation and the ChrY genetic element 
exerting regulatory properties. 

Studies on the Drosophila ChrY were the first to reveal its 
global gene regulatory properties and these studies have guided our 
investigation into whether a related mechanism was evolution- 
arily conserved by the mammalian ChrY. Research suggests that, in 
Drosophila, ChrY plays an important role in regulating the het- 
erochromatic state of autosomal and ChrX chromatin, thereby 
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Figure 5. Genetic variation in ChrY modifies the sexual dimorphism in 
myocarditis and leads to a sex bias in EAE severity. (A) Male and female mice 
from WT B6 and male mice from four B6-ChrY consomic strains were 
infected with 50 PFU CVB3 and their hearts were evaluated for myocarditis. 
Heart score of female B6 was subtracted from the heart score for male WT B6 
and each B6-ChrY consomic strain. Labeling on the x-axis represents the 
strain donating ChrY to B6. Results are representative of three individual 
experiments, n > 5 mice per strain. Significance determined by one-way 
ANOVA and Dunnett's multiple comparison test. (*) P < 0.05; (****) P < 
0.0001 . (B) Female B6 and male B6-ChrY consomic mice were immunized 
with MOG35_55 using the 2x protocol and the clinical score was monitored 
over 30 d. The significance of the differences in disease course among the 
strains was determined by two-way ANOVA (interaction [F= 1 .1 8; DFn = 31 9; 
DFd = 7770; P= 0.02], day post-injection [F= 1 32.8; DFn = 29; DFd = 7770; 
P< 0.0001 ], and strain [F= 23.5; DFn = 1 1 ; DFd = 7770; P< 0.0001 ]). (C) The 
significance of the differences in disease course between B6 female mice and 
each B6-ChrY consomic line was determined as described in Supplemental 
Figure S2 and grouped accordingly. The significance of the differences in 
sexual dimorphism among the groups was determined by comparing the 
best-fit values for the F < M and F > M grouping against the F = M group and 
by two-way ANOVA (interaction [F= 1 .1 4; DFn = 348; DFd = 8880; P= 0.04], 
day post-injection [F = 1 55.5; DFn = 29; DFd = 8880; P < 0.0001 ], and strain 
[f= 22.58; DFn = 12; DFd = 8880; P< 0.0001]) followed by Holm-Sidak 
corrected post-hoc multiple comparisons. (*) P < 0.05; (****) P < 0.0001 . 

regulating gene expression and piienotypic differences among 
males (Piergentili 2010). The degree of ChrY-linked regulatory var- 
iation inversely correlates with the number of ribosomal DNA 



(rDNA) units within the ChrY rDNA locus, which is a tandemly 
repeated array consisting of hundreds of units that comprise —10% 
on the entire ChrY (Paredes et al. 2011; Zhou et al. 2012). 

In the murine ChrY, >95% of the male-specific region is 
composed of hundreds of units of a tandemly repeated array 
(Alfoldi 2008). Each unit contains several gene copies from three 
testis-specific gene families, including Sly, Ssty, and serine-rich, 
secreted, Y-linked (Srsy) multicopy genes. The remaining 5% of the 
male-specific region is similar in structure and composition to the 
human ChrY and contains several single copy protein-coding 
genes, such as Sry, as well as another tandemly repeating array with 
dozens of units, with each unit containing a copy of the protein- 
coding gene Rbmy (Alfoldi 2008). While the function of Srsy and 
Ssty remains unknown. Sly represses the post-meiotic expression of 
ChrX and ChrY genes by maintaining repressive epigenetic marks 
in sex chromatin (Cocquet et al. 2009; Reynard et al. 2009; Ellis 
et al. 2011; Cocquet et al. 2012). Furthermore, Sly underlies the 
ChrY deletion mutation phenotypes in mice that exhibit a range of 
deficiencies in spermatogenesis depending on the extent of the 
deletion (Reynard et al. 2009). Rbmy is also important during 
spermatogenesis and encodes an RNA-binding protein that regu- 
lates mRNA splicing, and loss of RBMY expression may contribute 
to sperm abnormalities in mice and humans (Elliott et al. 2000; 
Szot et al. 2003; Skrisovska et al. 2007; Zeng et al. 2008; Navarro- 
Costa et al. 2010). 

We observed a highly significant correlation between the 
extent of copy number variation in Sly and Rbmy, but not Sstyl, 
with disease susceptibility. In addition, similar to the relationship 
between ChrY-linked regulatory variation with rDNA repeats in 
Drosophila, copy number variation in Sly and Rbmy is inversely 
correlated with the up-regulation of genes in immune cells. The 
regulatory properties of these ChrY multicopy genes are not likely 
exerted through their protein products since mRNA expression is 
testis and germ line specific (Elliott et al. 1996; Toure et al. 2005). 
However, we find our data to be consistent with a mechanism 
of chromatin remodeling as proposed for the Drosophila ChrY, 
whereby the multicopy genes may sequester proteins involved in 
chromatin dynamics (Dimitri and Pisano 1989; Lemos et al. 2008, 
2010). Therefore, a high frequency of the ChrY multicopy genes, as 
seen in ChrY^^ and ChrY^^°-^, may lead to the reduced availability 
of chromatin remodeling proteins, resulting in fewer euchromatic 
DNA regions and less transcriptional activity. Conversely, since 
ChrY^J^ possesses fewer copies of these multicopy genes, more 
chromatin remodeling proteins are freely available to exert their 
functions on chromatin, leading to more euchromatic regions and 
an increase in transcriptional activity. Furthermore, support for 
a mammalian ChrY-mediated chromatin modification mecha- 
nism stems from the surprisingly high number of alternative splice 
variants identified between the parental and ChrY consomic 
strains. Recent studies have identified a role for chromatin in the 
regulation of alternative splicing such that nucleosome position- 
ing along exons, plus the enrichment of chromatin-modifying 
proteins with splicing factor-binding capabilities, may regulate the 
recruitment of splicing factors to the pre-mRNA (Luco et al. 2011). 
Therefore, this study indicates that ChrY is a member of the reg- 
ulatory genome in males and supports an evolutionarily conserved 
mechanism of gene regulation by ChrY in the immune cells linked 
to Al disease. W^hether this is a direct consequence of copy number 
variation in these genes, or whether copy number variation is 
linked to other genetic features within ChrY, remains unknown. 
Consistent with our model, ChrY was recently linked to an NK 
and B-cell deficiency in mice, which is not due to changes in 
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protein-coding sequences as determined by exome sequencing, or 
due to expression differences in ChrY genes as determined by 
microarray and qPCR, but is rather due to a structural mutation 
arising from the deletion of one-third of the entire ChrY (Sun et al. 
2013). Therefore, this new study supports the structure of ChrY's 
long arm as the genetic element exerting gene regulatory prop- 
erties in immune cells. 

Nevertheless, additional studies must be conducted to effec- 
tively rule out the existence of genetic variation within the pseu- 
doautosomal region (PAR), an ~700-kb segment at the telomeric 
end of ChrY (Perry et al. 2001), that could also potentially in- 
fluence Al disease susceptibility in the ChrY consomic strains. In 
Mus musculus and Mus domesticus subspecies, the PAR contains one 
protein-coding gene, called steroid sulfatase (Sts), and the 3' end of 
the midline 1 gene (Midi), which is non-protein coding and spans 
the PAR boundary (White et al. 2012). Since the PAR is free to 
recombine with ChrX and exchange genetic information during 
male meiosis, all the ChrY consomic strains used in this study 
presumably possess the B6 PAR. However, the PAR boundary can 
regularly shift by several hundred kilobases due to the presence of 
repeat sequences present at the recombination site and, therefore, 
the PAR boundary may be different between the consomic strains 
(Kipling et al. 1996a,b). A dramatic example of a boundary shift 
arising in radiation-induced mutation rate studies at Oak Ridge 
National Laboratory is seen in ChrY^^^ ("Y-linked autoimmune 
accelerator") (Major 1954; Lane and Murphy 1972) whose PAR is 
extended by >4 megabases (Mb) and contains —19 additional genes 
that translocated from ChrX (Pisitkun et al. 2006; Subramanian 
et al. 2006; Santiago- Raber et al. 2008). However, we did not detect 
any significant expression level differences in genes residing within 
10 Mb of the PAR boundary in our microarray studies comparing 
B6-ChrY^J^ vs. B6 CD4^ T cells or SJL-ChrY^^^^ vs. SJL CD4-' T cells 
and macrophages, suggesting that the phenotypic differences ob- 
served in these consomic strains are not the result of a similar ex- 
pansion in the PAR. 

Only when we have obtained a detailed understanding of 
the regulatory genome will we be able to fully appreciate the 
mechanistic setting underlying complex disease traits. The recent 
unveiling of the ENCyclopedia Of DNA Elements (ENCODE) pro- 
ject has enlightened the scientific community on the vast wealth 
of information residing in non-protein coding DNA sequences 
(The ENCODE Project Consortium 2012). Our investigation into 
gene regulation by ChrY lends significance to what was once 
considered to be "junk" DNA and defies the belief that in mammals 
ChrY's contribution to sexual dimorphisms is limited to sexual 
differentiation, development, and associated male sex hormones. 
A comparison of the myocarditis heart score and the clinical EAE 
disease course of female mice with male B6-ChrY consomic strains 
indicates that sexual dimorphisms can be contextually dependent 
upon the cell type-specific genome-wide trans-eQTL properties of 
ChrY. Clearly, defining the epigenetic mechanisms whereby ChrY 
polymorphism influences the regulatory genome in human health 
and disease, and the role that such polymorphism plays in sexual 
dimorphism (Mank 2009; Parsch and EUegren 2013), is critical for 
the appropriate therapeutic management of complex diseases. 

Methods 

Mice 

All mice in this study were bred and maintained in the animal 
facility at the University of Vermont. The B6-ChrY consomic panel 



of mice used in this study was purchased from The Jackson Labo- 
ratory. The number of backcrosses preformed prior to cryopreser- 
vation is indicated next to each strain if this information was 
available on The Jackson Lab website. We then backcrossed each 
consomic strain to B6 at least one more time before conducting the 
experiments. C57BL/6J (B6), C57BL/6J-ChrY^2'^^^/^^^^J/NaJ (B6- 
ChrY^^''), C57BL/6JEi-ChrY^^^/^^J/EiJ (B6-ChrY^^^) N22, C57BL/ 
6JEi-ChrY^/J/EiJ (B6-ChrY^J) N17, C5 7BL/6JEi-ChrYL^^^^/EiJ (B6- 
ChrY^^^^^) N33, CSV^BL/G-ChiY^^^^^^/Vote] (B6-ChrYP^^) N14, 
C57BL/6JEi-ChrY^^/J/EiJ (B6-ChrY^^) N22, C5 7BL/6JEi-ChrY^JL/J/ 
EiJ (B6-ChrY^J^) N20, C57BL/6JEi-ChrY^™/EiJ (B6-ChrY^^) 
N20, C57BL/6JEi-ChrY^^^/^VEiJ (B6-ChrY^^^) Nil, B6Ei.MA- 
AChrY^^/^yJ/EiJ (B6-ChrY^^) N19, B6Ei.SWR-AChrY^™/J/EiJ (B6- 
ChrY^^) N20. SJL-ChrY^^^-^ mice were generated at the University 
of Vermont and backcrossed for a minimum of 15 generations and 
are available at The Jackson Laboratory (SJL/J-ChrY^^°/TeuJ). SJL/J 
mice were purchased from The Jackson Laboratory (Bar Harbor, ME). 
Animals were housed in specific pathogen-free conditions under 
National Institutes of Health guidelines, and all experiments per- 
formed in this study were approved by the Animal Care and Use 
Committee of the University of Vermont. 

Induction and evaluation of EAE 

Mice were immunized for the induction of EAE using the double- 
inoculation (2x) protocol and scored as previously described 
(Noubade et al. 2007) and for SJL (Spach et al. 2009). Histopatho- 
logical evaluations were done as previously described (Teuscher 
et al. 2004). 

induction and evaluation of myocarditis 

The H3 variant of CVB3 was made from an infectious cDNA clone 
as previously described (Knowlton et al. 1996). Mice were infected 
by intraperitoneal injection of 0.5 mL of PBS containing 50 PFU 
CVB3 and euthanized 7 d after infection. Hearts were fixed in 10% 
buffered formalin, paraffin-embedded, sectioned, and stained with 
haematoxylin and eosin. Sections were blindly evaluated on a 0-4 
scale, with 0 representing no inflammation, 1 indicating 1-10 le- 
sions, 2 indicating 11-20 lesions, 3 indicating 21-40 lesions, and 
4 indicating >40 lesions per section. 

Heart viral titers 

Hearts were homogenized in RPMI-1640 media (Mediatech) con- 
taining 2% ECS, penicillin, and streptomycin. Cellular debris was 
removed by centrifugation at 300^ for 10 min and supernatants 
subjected to a series of 10-fold serial dilutions in above media. Titers 
were determined by plaque forming assay as previously described 
(Knowlton et al. 1996). 

Serum testosterone 

Eight-week-old mice were bled by tail vein and the samples spun 
down for 20 min at 16,000^ to extract serum. Samples were stored 
at -80°C. Serum testosterone levels were measured using a testos- 
terone enzyme immunoassay (EIA) kit according to the manufac- 
turer's instructions (Assay Designs). Briefly, serum samples were 
diluted with 1 part steroid displacement reagent for every 99 parts 
sample. Each sample was diluted 1:10 and 1 : 20 in PBS and added to 
the 96-well plate in duplicate, followed by the primary antibody, 
and incubated at room temperature for 1 h while shaking. Wells 
were washed and conjugate was added to each well and incubated 
at room temperature for 1 h while shaking. Wells were washed and 
the pNpp substrate was added to each well and incubated at 37°C 
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for 1 h without shaking. Stop solution was added and the plate read 
immediately with optical density at 405 nm. 

FACS-pu rifled CD4^ T cells 

Three biological replicates were pooled from the axillary, brachial, 
and inguinal lymph nodes from five mice for each replicate. Cells 
were washed with PBS/1% PCS, pelleted at 300^ for 5 min at 4°C. A 
working stock of Live/Dead-775 dye (Invitrogen LlOl 19) was made 
by adding 1 fxL of reconstituted dye to 1 mL PBS and cells were 
resuspended in 100 ^xL of working stock per 1 X 10^ cells. Cells 
were incubated for 30 min on ice in the dark. Cells were washed 
with cold PBS/1% FCS, pelleted, and resuspended in 100 jxL of an 
antibody mixture containing 1 |jLg/mL PerCp-Cy5.5 anti-mouse 
TCRp chain and 1 ixg/mL Alexa Flour 488 anti-mouse CD4 anti- 
bodies in PBS/1% FCS. Cells were incubated on ice in the dark for 
30 min then washed with cold PBS/1% FCS, pelleted, and resus- 
pended at a concentration of 10 X 10^ cells per 1 mL of PBS/10% 
FCS. Cells were filtered through 50 fxm nylon mesh and sorted 
using BD FACSAria to collect ~3 X 10^ CD4"' T cells per sample. 

Adherent macrophages 

Adherent macrophages were isolated in culture from thioglycollate- 
elicited peritoneal exudate cells as previously described (Zhang 
et al. 2008). For each strain, three pooled biological replicates 
of adherent macrophages were obtained from five mice for each 
pool. 

RNA isolation 

RNA was isolated using the Qiagen RNeasy Plus Mini Kit following 
the manufacturer's instructions (Qiagen). RNA was eluted using 
30 |jlL RNase-free H2O and 0.5 |jlL Superase (Life Technologies) was 
added to each sample. RNA was quantitated on ND 1000 spectro- 
photometer v.3.3.1 (ThermoScientific) and quality assessed on 
a 2100 Bioanalyzer (Agilent). RNA was stored at -80°C and sub- 
mitted to the Vermont Genetics Network Microarray Facility. 

Mouse gene array target preparation 

Oligonucleotide microarray analysis of RNA expression levels was 
performed in the Vermont Genetics Network Microarray Facility 
using the Affymetrix GeneChip Platform (Affymetrix Inc.) ac- 
cording to manufacturer's protocols. In brief, the Nugen Ovation 
system v. 2 with SPIA RNA amplification was employed to convert 
50 ng of total RNA to cDNA. This isothermal RNA amplification 
system produces 5-12 fxg of anti-sense cDNA targets that is fol- 
lowed by several steps to produce sense strand cDNA to be frag- 
mented, biotinylated, and hybridized to the genechip. After pu- 
rification and fragmentation, biotinylated-cDNA targets were 
hybridized to the Mouse Gene 1.0 ST Arrays oligonucleotide arrays 
for 16 h at 45°C. Hybridized arrays were washed and stained with 
streptavidin-phycoerythrin followed by sequential incubations 
with biotin-coupled polyclonal anti-streptavidin antibody and 
streptavidin-phycoerythrin as a fluorescent amplification step. 
After staining, arrays were scanned (3000-7G Scanner, Affymetrix 
Inc.) and data collected for statistical analysis. 

Mouse miRNA array target preparation 

Total RNA containing small RNA was labeled using the FlashTag 
RNA labeling kit (Genisphere) and performed in the Vermont 
Genetics Network Microarray Facility according to the manufac- 
turer's instructions. Briefly, for each sample, 500 ng of total RNA was 



subjected to a tailing reaction (2.5 mM MnC12, ATP, Poly A 
Polymerase — incubation for 15 min at 37°C) followed by ligation 
of the biotinylated signal molecule to the target RNA sample (1 x 
FlashTag ligation mix biotin, T4 DNA ligase — incubation for 
30 min at 25°C) and addition of HSR stop solution. Each sample 
was hybridized to a GeneChip miRNA 2.0 Array (Affymetrix) at 
48°C and 60 rpm for 16 h then washed and stained on Fluidics 
Station 450 (Fluidics script FS450_0003) and finally scanned on 
a GeneChip Scanner 3000-7G (Affymetrix). 

Calculation of probe set statistics 

Raw GeneChip data (one DAT file for each chip) includes a col- 
lection of images, one for each probe and chip. Each image was 
summarized by Affymetrix GCOS software using one probe in- 
tensity (in CEL files, one per chip). Information from multiple 
probes was combined to obtain a single measure of expression for 
each probe set and sample. Probe-level intensities were calculated 
using the Robust Multichip Average (RMA) algorithm, including 
background correction, normalization (quantile), and summariza- 
tion (median polish), for each probe set and sample, as is imple- 
mented in Partek Genomic Suite, version 6.6 (Partek Inc., St. Louis, 
MO). Sample quality was assessed based on the relative log expres- 
sion (RLE), and normalized unsealed standard error (NUSE). Prin- 
cipal Component Analysis (PCA) was also used to look for outlier 
samples that would potentially introduce latent variation into 
the analysis of differential expression across sample groups. Linear 
modeling of sample groups and identification of alternatively 
spliced variants were performed using ANOVA as implemented in 
Partek Genomic Suite. The magnitude of the response (fold change 
calculated using the least square mean) and the P- value associated 
with each probe set and binary comparison were calculated, as well 
as step-up, adjusted P-value for the purpose of controlling the false 
discovery rate (Benjamini and Hochberg 1995). 

Bioinformatic analyses 

Data were analyzed through the use of the Mouse Genome In- 
formatics Gene Ontology Term Finder (www.informatics.jax.org/ 
gotools/MGI_Term_Finder.html) and by Ingenuity Pathways Anal- 
ysis (Ingenuity Systems, www.ingenuity.com). 

CD4^ T-cell activation and ELISA 

Three biological replicates of CD4^ T cells were purified from ax- 
illary, brachial, and inguinal lymph nodes pooled from five mice 
using the EasySep mouse CD4"^ T-cell enrichment kit and EasySep 
magnet following the manufacturer's instructions (Stemcell Tech- 
nologies). Nonpolarized effector T cells were generated as previously 
described (Noubade et al. 2007). Supematants were removed after 
24, 48, and 72 h and frozen at -80°C until cytokine levels were 
tested by ELISA. 

Quantitative PGR array 

Three biological replicates of CD4^ T cells were purified from ax- 
illary, brachial, and inguinal lymph nodes pooled from five mice 
using the EasySep mouse CD4"^ T-cell enrichment kit and EasySep 
magnet following the manufacturer's instructions (Stemcell Tech- 
nologies). High-quality RNA was isolated, tested for gDNA con- 
tamination, and hybridized to the RT^ Profiler 9 6- well PGR array 
(PAMM-086Z) run on Applied Biosystems 7900HT standard block 
real-time cycler by the Vermont Genetics Network DNA analysis 
facility at the University of Vermont. Data were analyzed using ABI 
Prism SDS 2.1 software (Life Technologies) followed by the RT^ 
Profiler PGR Array Data Analysis Version 3.5 (SABiosciences). 
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Gene copy number assay 

Tail DNA was isolated from a minimum of five male mice from 
each B6-ChrY consomic line and WT B6. DNA was quantified by 
Nanodrop and diluted to 5 ng/ixL in nuclease-free H2O. FAM dye- 
labeled PGR primers for Sly, Rbmy, and Sstyl were designed using 
the Custom Taqman Assay Design Tool on the Life Technologies 
website by entering the DNA sequence obtained from MGI. VIC 
dye-labeled mouse Tfrc Taqman Copy Number Reference Assay was 
used in the duplex real-time PCR reaction and run on the 7500 real- 
time PGR system according to the manufacturer's instructions. 
Data were analyzed using ABI Prism SDS 2.1 software followed 
by the GopyGaller Software Version 2.0 (Life Technologies). WT 
B6 was used as the reference strain and set to 10 copies for each 
of the genes. 

Data access 

The microarray data discussed in this manuscript have been de- 
posited in the NGBI Gene Expression Omnibus (GEO; http:// 
www.ncbi.nlm.nih.gov/geo/) (Edgar et al. 2002) under accession 
number GSE47440. 
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